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Three dimensional Ising model ferromagnets on different lattices with nearest neighbor interac- 
tions, and on simple cubic lattices with equivalent interactions out to further neighbors, are studied 
numerically. The susceptibility data for all these systems are analyzed using the critical Renormal- 
ization Group Theory formalism over the entire temperature range above Tc with an appropriate 
choice of scaling variable and scaling expressions. Representative experimental data on a metallic 
ferromagnet (Ni) and an elementary fluid (Xe) are interpreted in the same manner so as to estimate 
effective coordination numbers. 
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In the very extensive studies which have been devoted 
to critical phenomena attention has understandably been 
mainly concentrated on the regime in the immediate 
neighborhood of the critical temperature; data are an- 
alyzed using the leading terms in the Renormalization 
Group Theory (RGT) formalism. It is widely consid- 
ered that as correction terms proliferate outside this nar- 
row " critical region" they always lead ultimately towards 
Gaussian fixed point mean-field-like behavior above a 
temperature Tg determined by the Ginzburg criterion [ij . 
This criterion expresses a crossover from a fluctuation 
dominated critical regime to a high temperature Landau 
regime. Sophisticated theoretical and numerical studies 
have been made of this crossover 0-Q] in particular in the 
context of long range interactions. Experimental data on 
fluids, where effective interactions are expected to be long 
range, have been interpreted on this basis 

Here we discuss from a different perspective numeri- 
cal data on various 3d Ising systems, in particular mod- 
els where the interactions extend beyond nearest neigh- 
bor. We conclude that with appropriate observables 
and choice of scaling variable, for this family of mod- 
els at least there is no need to invoke a crossover or the 
Ginzburg criterion. The data can be convincingly inter- 
preted using the rigorous critical RGT formalism over the 
whole temperature range from to infinity. Applying 
Occam's razor, this approach is more economical concep- 
tually as it requires only the extended critical analysis 
but no separate analysis linked to a crossover. 

We discuss representative experimental data on a fer- 
romagnet and on a fluid, showing that they too can be an- 
alyzed in a transparent manner using the same approach 
and without invoking crossovers. We obtain quantita- 
tive estimates of the effective coordination number for 
the metallic ferromagnet Ni and for the elementary fluid 
Xe. 

Before RGT it was already established ll-l^ that 
in the fundamental scaling law for ferromagnetics the re- 
sponse parameter is the "reduced" susceptibility (see [H] 



for definitions) 



Xm - {{m {m)f) = Y,{S, ■ S,) ^ xt(/3)/xo(/3) (1) 



where the thermodynamic susceptibility 

Xt(/3) = [dm/dH]H^o 



(2) 



is normalized by the free spin susceptibility Xo(/3) cx j3. 
(As usual we will set the interaction strength J to 1 and 
write P = 1/kT). The critical behavior is written [isl - flij 



x(/3) = Txt(/3) « (1 + • • • ) 



(3) 



where e is an appropriately normalized scaling variable 
depending linearly on (T — T^) close to Tc- 

The thermodynamic ideal lattice gas analogue to xt 
is the isothermal compressibility on the critical isochore 
Kt ~ ~{dV / dp) /V so by strict analogy to the ferromag- 
netic case the parameter which should be scaled 13| is 
the compressibility normalized by the ideal gas compress- 
ibility {{N - {N)Y)/{N) where N is the total number of 
particles, i.e. 



TKt o^e-^[l + ■■■] 



(4) 



Because the gas-liquid order parameter is a scalar the 
fluid transition belongs to the Ising universality class 17 1 . 
Careful experimental measurements and analyses made 
over many years (see [lol |) have shown that the asymp- 
totic fluid critical exponents are indeed those of the short 
range Ising universality class. The real situation is how- 
ever more complicated than in the magnetic case because 
of the departure from vapor-liquid symmetry in real flu- 
thc fluid "susceptibility" is defined in 

In work based on the high temperature series expan- 
sion (HTSE) theory which was already firmly established 
in the 1950s, the critical scaling variable is taken to be 
either e = r [1 - or e = [1 - tanh(^)/tanh(^c)] 

11, 12, [l3l- Appropriate high temperature series are 
written rigorously as sums of terms where exact factors 
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multiply successive powers of e. However since the intro- 
duction of RGT, scaling expressions are often written in 
terms of e = i = [{T — Tc)/Tc]. This is just the simplest 
linear convention, but other scaling variables including 
T are just as legitimate as t in the region very close to 
Tc- (For instance in the special case of the square lat- 
tice Ising model an extremely sophisticated analysis uses 
(1/ sinh(2/3)-sinh(2/3))/2 as the scaling variable In 
addition r has obvious practical advantages in the tem- 
perature region well above Tc because its high tempera- 
ture limit is 1 (and not infinity as is the case for t). 

Systematic analyses of high temperature numerical 
data of nearest neighbor and long range interaction mod- 
els in the Ising universality class 0, 0, B ll^l , have been 
carried out with an effective temperature dependent sus- 
ceptibility exponent defined as 



7t(t) = -aiogxT(/3)/aiogt 



(5) 



following the phenomenological expression of Kouvel and 
Fisher |2l| . In the high temperature limit for any spin 
model XT — > /3 and i ^ T so 7((t) wiU automatically 
tend to 1 at high temperatures. There will necessar- 
ily be a crossover at an intermediate temperature where 
7t(i) passes from the critical 7 to 1. However, far from 
criticality the choice of variables is vital. Thus if the re- 
duced susceptibility parameter xW) rather than yt(/3) 
had been used for the definition of jt{t) in Rcf. [3| and 
following work, the effective exponent would have tended 
to at high temperature, not to 1. 

The full formal RGT Wcgner scaling expression for 
x(t) (rather than xt(0) the thermodynamic limit 
including confluent and analytic correction plus back- 
ground terms is written rigorously using r as [ij, Il6l , 



X{r) = Txt{t) = C^T-^[l + a^T'Fair) + b^rFtir) 



0-2- 



'-F^ir) 



(6) 



where 7, a, the confluent correction exponents 0^, and 
certain amplitude ratios are universal but the amplitudes 
themselves are not universal; the Fi are inflnitc analytic 
series in t normalized to 1 at r = 0. Because t — > 1 
as T — > 00 these developments remain well behaved at 
all temperatures above Tc, whereas because t diverges as 
T it is obviously very awkward to extend to high tem- 
peratures the analogous expression written in terms of 
t. 

The temperature dependent effective susceptibility ex- 
ponent deflncd in terms of r and x(r) in Refs. 1^ 24-27 
is 



1os{t) = -aiogx(T)/51og(T) 
with the equivalence 

7t = 7off(l - r) + r 



(7) 



(8) 



Very close to Tc the two effective 7 parameters are in- 
distinguishable but they have quite different properties 
as soon as the temperature difference increases. (Histor- 
ically it is of interest to note that the equation 



1/XT = {l/C)T{l - {Tc/T))y{Tc/T) 



(9) 



cited explicitly by Kouvel and Fisher |2l| to justify their 
analysis in the vicinity of Tc is precisely of the form of 
Eq. [51 with T as the scaling variable). 

At flrst sight the sets of infinite series of corrections 
in Eq. [6] appear rather forbidding. However, from in- 
spection of 5* = 1/2 high t emp erature series expansions 
(HTSE) (see for instance [l6|) there are exact closure 
rules at infinite temperature : C^(l -I- a,^ -|- ■ • ■ ) = 1 and 
7off(l) = z/3c 2J] where z is the coordination number 
and the (• • • ) represent the exact sum of all the higher 
order terms in Eq. [6] evaluated with r set equal to 1. 
It turns out that for the 3d Ising models which we will 
discuss explicitly, over a wide temperature region above 
Tc the leading [confluent] Wegner correction term domi- 
nates. As a convenient approximation all the remaining 
terms can be collected together into a single weak effec- 
tive correction term k^r^^ , giving a compact approxi- 
mate expression which can be used to fit the data over 
the entire temperature range above Tc : 



-) 



(10) 



so 



The closure rules then become 

C^{l + a^^k^)^\ 

and 

7 - {a^e + k^X^)C^ = z/3c 



(12) 



(13) 



Once the strictly defined critical amplitudes and 
are estimated for any particular model from data at tem- 
peratures close to criticality, k^ and are fixed also 
from the closure conditions, so the entire temperature 
dependencies of x(''') and of jesir) are determined. (An 
expansion of Eq. [TO] to include further explicit Wegner 
terms is possible when HTSE data and higher order cor- 
rection exponent values are available). 

While not rigorous except in the limits r ^ and 
T — i- 1, this ansatz as it stands already gives a representa- 
tion of the true behavior of x{t) which turns out to be ac- 
curate to the 10"'^ level over the whole temperature range 
above Tc for all the models we have studied. If 7oft(''") is 
transposed to 74 through Eq. [8| a crossover behavior re- 
sults (see Fig. 4); similar data have been analyzed using 
an approximant containing an implicit crossover function 

si 
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FIG. 1: (Color online) The normalized reduced susceptibil- 
ity xi''')'''"' a function of for the sc z26 (nearest, second 
nearest, and third nearest equivalent neighbors) Ising model. 
Lattice sizes L = 32, 24, 16, 12, 8 top to bottom (olive, pink, 
black, green, blue). The L > ^{L, (3) size independent enve- 
lope behavior region can be seen for each curve by inspection. 
The envelope fit curve, with extrapolation, is red 



The temperature dependent susceptibility x(/3,X) was 
evaluated on diamond, sc, bcc and fee lattices with near- 
est neighbor interactions, see (27l - [29j where the numerical 
techniques are described, and on sc lattices with equiva- 
lent interactions up to second, third, fourth or fifth neigh- 
bor following [25l . [30| . The coordination numbers are 
= 4, 6, 8, 12, 18, 26, 32 and 56 respectively. As all these 
models are in the 3d short range interaction Ising univer- 
sality class, in the analysis below the numerical values for 
the universal 3d Ising exponents were fixed at 7 = 1.2371, 



6 = 0.50 and v = 0.630 31, S^]. The critical inverse tem- 
peratures /3c for the various models were evaluated from 
the present data using the Binder cumulant g{l3, L) and 
the parameter VF(/3, L) introduced in |33|. The /3c values 
obtained from the finite size scaling analysis are in full 
agreement with previous estimates, in particular those of 
Ref. for the equivalent interaction models. The raw 
normalized reduced susceptibilities x{t^L)t'^ as a func- 
tion of for the 2:26 (nearest, second nearest, and third 
nearest equivalent neighbors) sc Ising model at different 
sizes L are shown as an example in Fig. 1. The envelope 
curve which can be seen by inspection corresponds to 
the thermodynamic limit (effectively infinite size) behav- 
ior. Data for each of the other models have qualitatively 
similar appearance (see for the nearest neighbor sc 
model). The susceptibility results for infinite L (extrap- 
olated for T close to zero) for the various coordination 
numbers z are exhibited in Fig. 2 in the form of plots 
of x(r)T''' against r^. The effective exponents 7eff(T) de- 
rived from these data are shown in Fig. 3. 



FIG. 2: (Color online) The normalized reduced susceptibil- 
ity xi''')''''^ ^ 8- function of for the diamond, sc, bcc, fee, 
zl8, 226, 2:32 and 266 sc Ising models, top to bottom (black, 
red, green, blue, cyan, pink, yellow, olive). The exponents are 
taken to be 7 = 1.2371 and 9 = 0.50. 



TABLE I: Values of the fitting parameters for x(''') in each of 
the models, as defined in Eq. 1101 
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0.3697 


1.245 


-0.1382 


-0.059 
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0.221655 


1.116 


-0.0914 


-0.0125 


3 
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0.15737 


1.038 


-0.073 


0.036 


1.6 


12 


0.102067 


1.022 


-0.062 


0.04 


1.5 


18 


0.06442 


0.871 


0.111 


0.038 


0.9 


26 


0.0430385 


0.765 


0.302 


0.009 


1 


32 


0.0343267 


0.703 


0.428 


-0.006 


3 


56 


0.0189291 


0.575 


0.800 


-0.06 


1.5 



The critical parameter estimates /3c (z), C^{z) and 
a^(z) and the approximate effective parameters k~^{z) 
and Ax(2:) are given in Table I. (The second correction 
term is always weak, so the values of A^(z) are not pre- 
cise as they depend very sensitively on the fit parameters 
chosen for the other variables). All the models, including 
those with longer range interactions, follow the critical 
scaling rules up to infinite temperature, with a gradual 
evolution of the critical amplitudes as z increases but 
without a trace of a crossover to mean-field behavior at 
high T . It is important to note that it is the coordination 
number z rather than the interaction range measured in 
terms of the nearest neighbor distance which is the key 
parameter (with weak lattice structure effects); the dia- 
mond, sc, bcc and fee arc all nearest neighbor lattices but 
they have significantly different values for C^^ and a^. 

There seems no obvious reason to expect a breakdown 
in these rules however large the range of interactions as 
long as there is a cut-off so that the range remains fi- 
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FIG. 3: (Color online) The effective exponent 7cff(T") as a 
function of for the diamond, sc, bcc, fee, zl8, z26, z32 and 
z56 sc Ising models, top to bottom (black, red, green, blue, 
cyan, pink, yellow, olive). 



FIG. 5: (Color online) The normalized reduced susceptibility 
Tx.t{t)t'' of Ni as a function of t* with the Heisenberg ex- 
ponents 7 = 1.349 and 9 — 0.55. The experimental data and 
units are taken from [36l (red circles) and [33] (black squares) . 
Following [33l and [39j. a small temperature independent term 
has been subtracted from the raw xt data. 




t/t 



FIG. 4: (Color online) The effective exponent 7t as a function 
of t/tx for the diamond, sc, bcc, fee, zl8,z26,z32 and 2:56 
sc Ising models, (black, red, green, blue, cyan, pink, yellow, 
olive), tx is defined following Ref. 



nitc; the correction amplitudes should continue to in- 
crease with increasing range. (If interactions fall off al- 
gebraically and sufficiently slowly, the models will leave 
the finite-range universality class [33. [35|). 

For comparison the data of Figure 3 translated appro- 
priately are shown (Fig. 4) in the form of a jt against 
t/tx where the normalization parameter is defined by 
tx = («x,z)~^^^ used in Ref. Q. It can be observed 
that there is no universality either in the position or the 



form of the individual curves. 

Experiments can be interpreted using the same ap- 
proach. The venerable experimental data for the sus- 
ceptibility of the ferromagnet Ni tabulated by Weiss and 
Forrer [36[ and by Fallot 37 1 are exhibited in Fig. 5 in 
the same form as that used for the numerical data in Fig. 
2. Here we consider Ni as a 5 = 1/2 Heisenberg local mo- 
ment system and so use the Heisenberg exponent values 
7 = 1.396 and = 0.55 Following Fallot himself 

and [1^, we have subtracted out a small temper- 



37 



ature independent susceptibility term, which could well 
come from an orbital contribution (see [i^l for the case 
of Co). As in the Ising models with higher coordina- 
tion numbers shown in Fig. 2, the normalized reduced 
susceptibility increases almost linearly with over the 
wide range of temperatures covered which extends to 3Tc, 
i.e. ^ 0.8. The ratio between the asymptotic critical 
value of x(t)t'^ and the estimated extrapolated infinite 
temperature value (equal to 1 for spin 1/2 in the appro- 
priate units) can be taken as a measure of the effective 
C-^. For Ni the observed ratio is about 0.60, or alterna- 
tively the correction amplitude 0.65. 

As an example of a gas-liquid transition we consider 
the susceptibility xt (defined as the derivative of the 
density by the chemical potential xt = {dp/dij)T) of 
Xe on the critical isochore for the liquid-gas transition, 
for which results from careful experiments based on light 
scattering techniques are tabulated in Ref. 
6 these data are plotted in the form Txtt'' 
with the 3d Ising exponents. (It can be noted that the 
susceptibility xnn defined by Orkoulas et al [l9j in their 



41[. In Fig. 

against r^, 
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FIG. 6: (Color online) The normalized reduced compressibil- 
ity TKtt'^ of Xe as a function of r^. Experimental data and 
units taken from [ZH ]. 



FIG. 7: (Color online) The effective critical correction ampli- 
tudes ax{z) for the Ising models as a function of coordination 
number z, Table I, plus points for Ni and Xe obtained from 
the plots of the experimental data in Figs 5 and 6 (arrows). 



analysis of the hard-core square- well fluid is Txt] Ork- 
oulas et al also use t as the sealing variable). Again 
the figure shows an essentially linear increase of Txtt'^ 
with just as in the numerical plots for the Ising mod- 
els with large coordination numbers. The ratio of the 
critical limit to the extrapolated high temperature limit 
is Cy^ ~ 0.55, or alternatively the correction amplitude 
~ 0.75. This value is broadly consistent with the val- 
ues a^^ = 1.3(2) [4l| and = 1.08 [1] estimated from 
previous analyses based on the same data set but using 
different scaling rules. The complications associated with 
the asymmetry in the fluid phase diagram should be kept 
in mind, but this plot suggests that as in the magnetic 
case even if fluid data were available to much higher r 
within the present approach there would be no need to 
invoke a crossover to mean-field like behavior. In Fig. 
7 the values of for the numerical models are plotted 
against z, and the effective values for Ni and for Xe are 
indicated by arrows. From this figure we can estimate 
the effective coordination number z for Ni and for Xe. 
To obtain a more quantitative estimate the Ni experi- 
mental data should be compared to numerical results for 
Heisenberg spins on an fee lattice for different z rather 
than for Ising spins on an sc lattice. From data compar- 
ing 7(t) on Heisenberg spins on sc, bee and fee lattices 
2J] it appears that the numerical a-^{z) plot will be of 
similar form for Heisenberg spins as for Ising spins. One 
can then estimate that approximately z ~ 45, or in other 
words the effective interactions extend to between two 
and three lattice spacings. Obviously for real physical 
systems the model of equivalent interactions with cut- 
off is only a rough approximation to the true situation 
but the effective z is a useful indicative phenomenological 



parameter. 

For Xe the effective corresponds to an effective se 
Ising coordination number z ^ 55. When data for some 
other fluids are expressed graphically in terms of 7e/ / (t) 
Q, 7e// initially increases slightly with increasing tem- 
perature, meaning that is negative. (The experimen- 
tal fluid data sets generally extend only over a narrow 
temperature range above Tc] for small r, jt(t) is practi- 
cally indistinguishable from ^rir) and values estimated 
for are not sensitive to which scaling expression is 
used) . Negative a^, values have been observed in all aque- 
ous electrolyte solutions and also in many non-aqueous 
ionic solutions It has been suggested [1] that as a 

general rule simple fluids have positive and complex 
fluids negative a^. A comparison with Table I and Fig. 
3 indicates that for the negative systems the effective 
coordination number z is 12 or less while the positive 
systems have much higher effective z values. 

The effective coordination number is a fundamental 
parameter for understanding the magnetism of metallic 
ferromagnets which can often be considered either from a 
band or from a local moment perspective. In the case of 
liquids, it should be possible to make a systematic clas- 
siflcation of effective coordination numbers and to link 
these z values to the interatomic potentials used for cal- 
culating structure functions. 

In conclusion, the analysis given above leads to a sim- 
ple overall physical scenario in which for a family of 3d 
Ising models the temperature dependence of the reduced 
susceptibility over the entire temperature range from Tc 
right up to inflnite temperature is explained using the 
critical RGT formalism with appropriate Wegner correc- 
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tions and without the need to invoke a restricted "crit- 
ical region" or any form of high temperature crossover. 
The approach is conceptually economical and leads to 
a transparent interpretation of the differences in behav- 
ior from model to model and from system to system 
within a universality class; there is a strong correlation 
between the coordination number and the strength of the 
non-universal amplitude of the leading confluent Wegner 
term, which dominates the corrections. 

We acknowledge gratefully an interesting discussion of 
the experimental high temperature effective exponents 
with Ralph Chambcrlin, and helpful explanations on the 
thermodynamics of fluids from Jan Sengers. This re- 
search was conducted using the resources of High Perfor- 
mance Computing Center North (HPC2N). 
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